Explanation for the acceleration of the universe's expansion without dark energy 
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We have observed the acceleration of the expansion of the universe. To explain this phenomenon, 
we usually introduce the dark energy (DE) which has a negative pressure or we need to modify 
the Einstein's equation to produce a term which is equivalent to the dark energy. Are there other 
possibilities? Combining our previous works of statistical mechanics of self-gravitating system with 
the derivation of van der waals equation, we propose a different matter's equation of state (EoS) in 
this paper. Then we find that if the matter's density is low enough, its pressure can be negative, 
which means that it is the matter that drives the expansion's acceleration. So here we will not 
need to add the DE to the universe. Our results also predict that the universe finally tends to be 
dominated by an approximate constant energy density, but its value can be smaller than DE. The 
data of Supernova can not differentiate our model from the standard model, but they may indicate 
some deviations from AC DM. 

PACS numbers: 



I. INTRODUCTION 

Under the assumption of isotropy and homogeneity 
of our Universe at large scale structure, people made a 
successful explanation to the cosmology accelerating ex- 
pansion [1] by introducing a new component called DE, 
which constitutes of about 70% energy of the universe. 
Up to now, the most successful view about the DE is 
the AC DM model, in which the energy density of DE is 
constant. Although ACDM model is consistent very well 
with all observational data, it faces the fine tuning prob- 
lem [2]. Alternatively, plenty of other DE models have 
also proposed [3-6] , but almost all of them solve the ac- 
celeration problem either by introducing new degree of 
freedom or by modifying gravity, which are challenges 
both in Cosmology and in Nuclear Physics. 

However, do we really need the DE? Or do we need so 
much DE? We notice that commonly the matter's EoS 
comes from the classical thermodynamics without consid- 
eration of the gravity, has it been confirmed to be reason- 
able at the cosmological scale? We think that the answer 
may be not sure, and our latest studies have provided 
different understandings of the long-range statistical me- 
chanics. In fact, if we apply the Boltzmann entropy 



S = - /ln/dT 



(1.1) 



into the self-gravitating system, we will always obtain the 
infinite mass and energy with the principle of maximum 
entropy [7], which is a serious problem. [8] preliminarily 
studies the entropy form taken by [9] and proposes that 
a self-gravitating system's entropy is at a saddle point 
and not a maximum. [10] completes the variation pro- 
cess of entropy and confirms that we can obtain an EoS 
which is different from the classical one but can solve the 
problem we mentioned above. In [11] we compare this 



EoS with the van der waals equation and speculate that 
our EoS should be a requirement of long-range statistical 
mechanics, which is also shown in this paper. 

In this paper, we think that the general relativity can 
also have similar effects on the matter's EoS at the cos- 
mological scale. So in the next section we will show our 
findings and propose a reason for the universe's accelerat- 
ing expansion, which will not need the DE. In section III, 
we will use the data of supernova to shorten the range of 
the values of the parameters in our model. Some discus- 
sions and conclusions will be made in the final section. In 
this paper we set c = 1 and use "0" to denote quantity's 
current value except Q x o below. 



II. EXPLANATION FOR THE ACCELERATION 

According to the cosmological principle, the universe is 
homogeneous and isotropic. Commonly the matter can 
be treated as the ideal fluid whose energy-momentum 
tensor can be written as 



= {p + p)U»U" + pg p 



(2.1) 



where p is the energy density, p is the pressure, and C/ M is 
the velocity four-vector. When we substitute it into the 
Einstein's equation, we can get the Friedmann equation: 

8ttG , 



AnG(- + p)a,a 2 + n 



-pa 



(2.2) 



where a is the scale factor. The energy-momentum con- 
servation law is 



3d . , 

P = (P + P), 

a 



(2.3) 



which also can be obtained from eq.(2.2). Combining the 
EoS p = p(p) with eq.(2.2), we can get the solution of a, 
p and p. After the inflation, the universe is assumed to 
be a gas whose pressure is [12] 
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-p. 



(2.4) 
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where v 2 is the mean-square velocity, m is the mass of 
particle, and T is the defined temperature. We usually 
use p — and p = p/3 to denote the matter (v 2 — 
« 1) and radiation's (v — 1) EoS respectively, but 
by eq.(2.4) we can not understand the nature of the DE 
whose pressure is negative, although we commonly treat 
that the DE does negative work to the universe. 

Notice that eq.(2.4) can be derived by two ways [13]: 
we can use the formula of the classical statistical mechan- 
ics to calculate the matter and radiation's EoSs, both of 
which satisfy eq.(2.4); it also can be obtained by the ki- 
netics alone, which needs to consider the bombardment 
by the particles (rigid body) on the walls of the container 
to calculate the kinetic pressure. However, the classical 
statistical mechanics does not consider the effects of grav- 
ity. Besides, the wall of the container does not exist in the 
universe, and in fact we can treat that the wall is equiva- 
lent to a potential well with minus infinity at the bound- 
ary for the free particles, while the gas in the universe is 
in the gravitational potential well which is different from 
the wall. If we treat the universe as a thermodynamical 
system where the particles interact by gravity, will these 
two ways mentioned above still be available under the 
cosmological background where the effect of general rel- 
ativity can not be neglected? We think that the answer 
is yes for the radiation, because the effect of gravity can 
be neglected compared with the electromagnetic force; 
but for the matter whose main component is the dark 
matter where the gravity plays the most important role, 
the answer may be not sure. We consider this problem 
by two sides. On the one hand, in [10] by the method of 
statistical mechanics we have obtained an EoS which is 
used to describe the thermodynamical equilibrium state 
of a virialized self-gravitating system: 



p = /3p + ap 4/5 , 



(2.5) 



where (3 and a are positive, and (3 has been identified 
to m/kT in [11]. Notice that here we take the Newto- 
nian gravity and the power index of the second term of 
eq.(2.5)'s right ride must be less than 1 to ensure of the 
finite mass and energy of the system. From eq.(2.5) we 
can obtain p = p(p), which can be approximately written 
as [11] 



kTp ,kT .9 4 

p = a{ — ) 5 p 5 . 

m m 



(2.6) 



Of course, eq.(2.6) is different from eq.(2.4) even we let 
v 2 = 3kT/m, and some different results possibly can be 
obtained. On the other hand, we consider the van der 
waals equation: 



n 2 a N 



(p+y^)(v-nb) 



NkT, 



(2.7) 



where a and b is used to describe the attractive and re- 
pulsive interactions among the molecules of gas respec- 
tively, and they are determined by the potential of the 



molecules' interactions by the method of cluster expan- 
sion in the canonical ensemble [see some books of sta- 
tistical mechanics such as 13], so we emphasize that the 
existence of a and b is determined by the statistical me- 
chanics. If the potential of the molecules' interactions 
has other forms, a and b may still appear in eq.(2.7)but 
their values and other places may be changed. Especially 
when we set V = Nm/p and b = 0, eq.(2.7) becomes 



kT , n 

P= — P ~ a \~jTj 
m Nm 



)V, 



(2- 



and is very analogous to eq.(2.6), which makes us specu- 
late that existence of a in eq.(2.6) also is a requirement 
of possible statistical mechanics which can deal with sys- 
tems with long-range interaction of the Newtonian grav- 
ity, and the change of index from 2 to 4/5 may be caused 
by the change from the interactions among molecules to 
the Newtonian gravity. Then what is the case of that the 
particles are in the gravitational field described by gen- 
eral relativity? The answer is that such theory still does 
not set up now, but we may do some guess of its results 
as the following: based on above analysis, we find that 
/? is related to the temperature, a is necessary to denote 
the effect of the long-range statistical mechanics, and nei- 
ther of them will disappear if the gravity changes from 
the Newtonian to the general relativity; but the index 
4/5 in eq.(2.6) may be changed in the general relativity. 
So we assume that at the cosmological scale the EoS of 
the matter may be written as 



Pn 



-Pn 



where po = §^y, t is a constant, and v and £l x o are about 
constants or very slowly changing functions of time at 
the matter dominated era. Here we temporally do not 
assume that v 2 is as small as its value in KCDM . Notice 
that the form of eq.(2.9) can be directly borrowed from 
the van der waals equation, and only t < 1 is speculated 
from our previous works. Eq.(2.9) can be named "the 
statistical matter model" . 

Now we will analyze the effect of eq.(2.9) on the uni- 
verse's evolution. Our model only modifies the pres- 
sure of matter, so the cosmic history before the matter- 
dominated era is believed to be littlely changed and here 
we only consider the later time. From observations we 
know that the unverse is always expanding, so p m is de- 
creasing with time. When p m is large enough, p m is pos- 
itive but is smaller than eq.(2.4). Besides, there exists an 
era when p m ~ like KCDM. But it is evident that the 
pressure in eq.(2.9) can be negative If 



(2.9) 



Po v 2 



(2.10) 



which tells us that when the matter reaches the thermo- 
dynamical equilibrium, its pressure will be negative if its 
density is low enough. This really is a surprising thing 
and need to be further confirmed, but it may be natural 



3 



in our model. If the matter's density is always decreasing, 
it will always can satisfy eq.(2.10), then can this negative 
pressure cause the Unverse's accelerating its expansion? 
From eq.(2.2) this require that p < —p/3, i.e. 



Po 



3tt x Q 



3v 2 



1 



(2.11) 



Of course this equation will also can be satisfied because 
of the expansion of the universe, so we find that if the 
long-range statistical mechanics is really different from 
the classical one and can produce an EoS like eq.(2.9), it 
will be the matter that drives the universe's accelerated 
expansion! The next question is that how the universe 
will evolve in the future, which requires some detailed cal- 
culations. We substitute eq.(2.9) into the eq.(2.3), then 
the matter's density will evolve as 

Prn = { C{fr B ^+n xQ ^ 
PO 



B 



- [(|) i - t (^r+-+(^)^] (2.i2) 



where 



B = 1 + — , C = B(l - tiro) 1 - 1 - O x0 , (2.13) 

and Sl r is the current radiation's density parameter. No- 
tice that 1 < B < 4/3. From observations we know k ~ 
and f2 r0 ~ 0.04, if we define 



= — 

PO 



(2.14) 



its current value f2 m o = 1 — ^ro will be much closer to 1, 
which is a major difference of our model from ACDM . 
When Q x0 = (or a is small enough so that Q x o can 
be neglected) and B tends to be 1, eq.(2.12) will become 
p oc a~ 3 . If B = 1 and t = 0, the form of cq.(2.12) will 
be the same as a direct addition of the matter and DE, 
so our model includes the result of ACDM, but we still 
provide a different understanding of the universe and any 
deviations from B = 1 or t = will indicate the prefer- 
ences of our model. Because of the universe's expansion, 
from eq.(2.12) we can easily find that the density will 
finally approximately become a constant 



Px = Po(-g-) 1 -* 



(2.15) 



If we neglect the variations of v and Sl x o with time. So 
our model predicts almost the same destiny of the uni- 
verse as ACDM: the standard model says that the en- 
ergy consists of the radiation, the matter and the DE, 
with the expansion of the universe the DE will dominate 
and the expansion will accelerate; while our result states 
that because there is some unknown physics such as long- 
range statistical mechanics, there is a no-zero minimum 
of the energy density of the matter even the universe will 
always expand, which also accelerates the expansion of 



the universe. However, this minimized density can be 
smaller than the DE, because: from cq.(2.2) we know 
that a > requires that p decreases slower than a -2 , so 
from eq.(2.12) we know that not only the constant term 
but also other terms whose power indexes are larger than 
-2 make the expansion accelerate. So compared with DE, 
the value of eq.(2.15) should be smaller if t > 0. Then if 
in the future the universe's density becomes a constant, 
according to eq.(2.2) our results may give a smaller value 
of the Hubble parameter than the standard model, which 
also differentiates our results from the standard model. 

Before analyzing the data of SNIa, we need to calcu- 
late the Hubble parameter. It can be expressed by four 
parameters (t, Q x o, B, ho) in our model, 

H(z) _ ,c(i + z)- 3B ^ + n x0 .^_ 

- { 5 J 1_t 



H y B 

+n r0 (i + zf. 



(2.16) 



where ho is the present value of the Hubble parameter 
H in unit lOOkm/s/Mpc and z is the redshift. 



III. COMPARED WITH THE DATA OF SNIA 
A. Data Analysis 

To check our proposals made in Sec. II and to con- 
strain the model parameters, we take use of the Union2 
dataset [14] and the Hubble evolution data. 

The Union2 dataset contains 557 type la SN data and 
covers the redshift range z — [0.015, 1.4], including sam- 
ples from other surveys, such as CfA3 [15], SDSS-II Su- 
pernova Search [16] and high-z Hubble Space Telescope. 

We fit the SNIa data by minimizing the % 2 value of the 
distance modulus. % 2 „ for SNIa is obtained by compar- 
ing theoretical distance modulus pth(z) = 51og 10 [(l + 
z) J z dx/E(x)} + Mo (po = 42.384 - 5 log 10 ho, E{z) = 
H(z)/H z ) with observed /i b of supernova: 



557 



Xsn ^ ] 



[Pth{Zi) ~ Pob{Zi)f 



To reduce the effect of no, we expand xtn with respect 
to (io [17]: 



X 2 sn = F + 2G(i Q + H(i 2 (3.1) 



where 



F = E 

i 

G = J2 



[pth(zj;po = 0) ~ Pob(zj)]' 2 
o\zi) 

Pth{zi,po = 0) - (iob(zi) 



<J 2 {Zi) 



4 



t 




B 




v-AOl +0.346, +0.562 


n «9/i~ u ' luv - - ( ^ lbl 

U.O^1 +0 200, +0.271 


1.030 


„ _„ . -0.02b, -0.U52 
U - ' /4 +0.050, +0.074 


0.051 


0.734 


1.008 


0.738 



TABLE I: Expected values, 1 — a and 2 — a error of 
t, Quo, B, ho in this model. The second line is in form of 
Expectation"}^' ^2ct- And the best-fit values are also shown 
in the last line, which are different from the expected values. 



(3.1) has a minimum as 



X 



sn.mtn 



F - G 2 /H 



which is independent of /iq. In fact, it is equivalent to 
performing an uniform marginalization over /io, the dif- 
ference between x 2 sn an d the marginalized xin i s J us t a 
constant [17]. We will adopt xtn as the goodness of fit 
between theoretical model and SNIa data. 

We also use the 12 Hubble evolution data from [18] 
and [19], its Xh ^ s defined as 



2 _ ^ [H(zj) - H ob { Zi )f 
Xh - 2__ 

i=l 



07 



Note that the redshift of these data falls in the region 
z e (0,1.75). 
In summary, 



2 _ ~2 2 
Xtot XsnXH 



and we assume uniform priors on all the parameters. We 
also assumed the prior that the age of our universe To 
satisfies lOGyr < T < 20Gyr. 




1.00 1.0S 1.10 1.15 1.20 1.25 1.30 




FIG. 1: ID marginalized distribution probability of 
t, ft x o, B, ho. 
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B. Results and Discussions 



FIG. 2: 68% and 95% contour plot in t - B plane. Red dot 
in the center is the expected value. 



The analysis is performed by using the Monte Carlo 
Markov Chain in the multidimensional parameter space 
to derive the likelihood. Naturally, we employ some phys- 
ically obvious limitations to make the estimation of pa- 
rameters more robust, for example, we set t < 1, 1 < B < 
4/3 in our analysis. We first investigate the constraint on 
the model parameters, and the best-fit values and errors 
of parameters are summary in Tab. I. We also plot the 
ID marginalized distribution probability of each param- 
eter, shown in Fig.l. The likelihood distribution shows 
a remarkable deviation from the Gaussian distribution, 
which result in a discrepancy between the best-fit values 
and the expected values of the parameters, as consistent 
with the result in Tab. I. 

Further, we focus on parameters {t,B), which play a 
more important role in our model. The 2D contour plot is 
shown in Fig. 2, which shows that the results of standard 
ACDM model is contained in our model. The 2a errors 
mean that the case with t > 0.80 is excluded at 95.4% 
confidence level, which is consistent with our assumption. 
But constraint on B can not be obtained. 



From above results, we find that the data of SNIa 
can not differentiate our model from ACDM. We need 
to consider the perturbation theory to further test our 
model, which will be another work. But from Tabled 
and Figure. 1 we can see that the data seem to indicate 
some deviations from B — 1 and t = 0, i.e. ACDM. 
Notice that according to current criteria, the dark mat- 
ter is called to be cold if v < 0.1 , which means that 
B < 1.0033, so there is only a very short range of B for 
the cold matter, however, both the best-fit and expected 
value of B are larger than 1.0033, which suggests that 
the matter may be not so cold as ACDM predicts. This 
conclusion also agrees with some works [20, 21] which 
aim to use the warm dark matter to solve the contradic- 
tions between observations and numerical simulations of 
ACDM at the galaxy scale. 



1 see http://en.wikipedia.org/wiki/Dark_matter 
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IV. CONCLUSION 

Our previous works have proved that the thermody- 
namics of gravitating systems at the galaxy scale may be 
different from the classical one, and we obtain an EoS like 
eq.(2.6) which is in analogy with the van der waals equa- 
tion. So we speculate that it is the long-range statistical 
mechanics that produces the eq.(2.6). Based on possible 
effect of the general relativity on the matter's EoS, we 
propose that at the cosmological scale the form of the 
matter's EoS may be like the eq.(2.9). Then we find that 
with the universe's expansion the density will decrease 
with time to be a no-zero minimized constant matter 
density, and it is the matter that drives the expansion's 
acceleration. So our model provides a new kind of ex- 
planation for the expansion's acceleration, but there are 
three major different results of our model from ACDM: 

(1). TO = 1 — fl r is close to 1 and is much larger than 
its value in ACDM; 



(2) . B and t may be not exactly to 1 and respectively; 

(3) . In the future when the unverse's density becomes 
a constant, its value can be smaller than DE. 

It is more easy to compare our model with ACDM 
by the second point, so we use the data of SNIa to con- 
strain the value of B and t. From the results we can not 
differentiate our model from ACDM, but there seems 
some indications of deviations from ACDM, especially 
we find that the matter may be not so cold as previously 
believed. To further test our model, we need to consider 
the density's perturbation in the future. 
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